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ABSTRACT 

Aims. The distribution of angular momentum of the matter during X-ray bursts on neutron stars is studied by means of 3D axi- 
symmetric general relativistic hydrodynamics. 

Methods. The set of fully general relativistic Navier-Stokes equations is solved implicitly using the implicit solver GR-I-RMHD in 
combination with a third order spatial and second order temporal advection scheme. The viscous operators are formulated using a 
Kerr-like metric in the fixed background of a slowly rotating neutron star whose radius coincides with the corresponding last stable 
orbit. The importance of these operators and their possible simplifications are discussed as well. To verify the consistency and accuracy 
of the solution procedure, the time-dependent evolutions of non-rotating heat bubbles during their rise to the surface of a white dwarf 
are followed and compared with previous calculations. 

Results. In the rotating case and depending on the viscosity parameter, o-mr, it is found that the viscously-initiated fronts at the center 
of bursts propagate at much faster speed than the fluid motion. These fast fronts act to decouple angular momentum from matter: 
angular momentum is transported outwards while matter sinks inwards into the deep gravitational well of the neutron star, thereby 
enhancing the compression of matter necessary for initiating ignition, that subsequently spreads over the whole surface of the neutron 
star on the viscous time scale. Based on the numerical simulations, we find that a viscosity parameter o-m,. = (9(0.1) is most suitable 
for fitting observations of neutron stars during X-ray bursts. It is argued that the spin up observed in the cooling tails of X-ray bursts is 
a transient phase, which eventually should be followed by a spin down phase. This delay can be attributed to a significant lengthening 
of the viscous time scale due to rapid cooling of matter in the outer layers. 
Conclusions. 

Key words. General relativity: neutron stars -black holes - X-ray bursts. Methods: numerical - hydrodynamics - relativistic 
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^ ' 1. Introduction 



Part of the observed neutron stars (NSs) belong to the family 
of ultra-compa ct objects, in which genera l relativistic effects 
are prominen t (IShapiro & TeukolskyTUgsl IStergioulasl l200l 
|Psaltisll2008[) . Depending on the equation of state, NSs may live 
even inside their last stable orbits, making the conversion ef- 
ficiency of gravitational energy into radiation eve n larger than 
that o f accreting Schwarzschild black holes (BHs, ICamenzindl 
120071) . In low mass X-ray binaries, relativistic jets have been 
also observed to emanate from around accreting NSs with bulk 
Lorentz factors tha t are comparable to those emanating from 
around black holes (lMigliarill2006ll2008h . 
On the simulation site, numerous general relativistic hydro- 
dynamical calculations have been carried out to study dif- 
ferent aspects of NSs, such as format ion, merger , inner 
structure, ac cretion or j ets around NSs ('T hielemannl I199C : 
Liebendorfer et alll2002l:lM arti & Muller, 2003 1; lOzel & Psaltis . 
2003; Shibata', '2003; 'McKinnev, 2006; Shibata etaLll2006 ; 
Abdikamalov et al., 2008; Zachariah et al., .2008. see also the 
references therein) 
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iDuez et"aL 1 (12004 carried out general relativistic calculations to 
study the formation of hypermassive NSs, taking into account 
the effect of viscosity. The authors found that viscosity drives 
the NS's inner core into rigid rotation and simultaneously trans- 
ports angular momentum outwards into the outer envelope. As a 
consequence, the core is found to contract in a quasi-stationary 
manner while the outer layers expand to form a differentially ro- 
tating torus. 

This behavior is similar to accretion of matter via rotating disks. 
Here the viscosity acts to decouple matter from angular momen- 
tum, in that it transports angular momentum outwards, while 
forcing the matter to sink deep er into the gravitational well of 
the central obiect( ftinglellT98l1) . In the absence of viscosity, an- 
gular momentum as well as magnetic fields in ideal MHD are 
frozen-in to the matter. Thus, while strong magnetic fields are 
essential to enable rigid rotation, viscosity, on the other hand, 
drives the outer layers into differential rotation. 
The type of rotation in the outer layers, differential or rigid rota- 
tion, may have profound effects on the conditions leading to the 
X-ray bursts observed on NSs jBildsten & Strohmaveil [19991; 
ISpitkovskv eTail l2002h . 

Indeed, recent observations of X-ray bursts revealed the so 
called burst oscillations, in which a spin-up or spin-down of 
the NSs in their cooling tails have been detected, reac hing a 
plateau on the asymptoti c hmit (see Strohmayer etdH 119971 ; 
IStrohmaver & Markwardi 119991; [Sttohmaver. ,20011 for a de- 
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tailed discussion). It has been argued that the increase/decrease 
of the spin of NSs during bursts is connected to the redistribution 
of ang ular momentum of the thermonuclear shell (Strohmaver, 
1200 Ibl) . Accordingly, when a thermonuclear shell starts to ex- 
pand at the burst onset, the moment of inertia increases while 
its spin decreases. When the shell starts to contract and subse- 
quently recouples to the NS, the inertia decreases and the spin 
increases. 

Also, several X-ray bursts on NSs display a spin do\y n rather 
than spin up in their cooling tails ('StrohmayeJ [19993). In this 
case, however, it was suggested that the spin down probably be- 
gins anew episode of thermonuclear energy release, most likely 
in layers underlying those re sponsible for the initial runaway. 
ICumming & BildstenI (l2003b . however, investigated in detail the 
hydrostatic expansion during bursts and the expected change of 
spin due to angular momentum conservation and concluded, that 
a shell expansion/contraction alone cannot explain the mecha- 
nisms underlying the observed spin-up/down of the NSs during 
bursts. Noteworthy is that the model in which ignition starts at 
a point and spreads over the whole surface of the NS via burn- 
ing fronts appears to fit observations, which reveals that the X- 
ray e mitting area increases during the bursts ( Strohmaver et all 
Il997i) . However, the role of rotation, the nature of these burning 
fronts and the manner they affect their surrounding are poorly 
understood. 

In this paper we present a first attempt to model the rotational 
evolution of thermally induced bursts beneath the atmosphere of 
a rotating NS and to study the role and effects of the viscosity 
on the redistribution of angular momentum under strong gravi- 
tational field conditions. Our investigation relies on employing 
a general relativistic hydrodynamical solver, in which turbulent- 
eddies have the effect of friction that gives rise to an enhanced 
re-distribution of angular momentum. 

The paper runs as follows: In Section 2 we describe the ad- 
ditional viscous operators that have been incorporated into the 
solver to study the viscous-redistribution of angular momentum. 
The results of several model calculations aimed at studying the 
distribution of angular momentum during X-ray bursts on NSs 
are presented and discussed in Section 3, while in Sec. 4 the re- 
sults are summarized. 



2. The general relativistic Navier-Stokes equations 

The set of the general relativistic hydrodynamical equations and 
their derivations are well described in Sec. 2 of Hujeirat et al. 
(2008). In this section we list the viscous operators of the mo- 
mentum equations, which we have incorporated into the implicit 
solver. 

The stress energy tensor for viscou s flows has the following f orm 
jRichardson & Chund. l2002t iFond [20031: ICamenzindi l2007l) : 
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where fi, v are indices that correspond to the four coordinates 



{t, r, 6, ip} and Tpp, T^j^ denote the stress energy tensor due to 
perfect and viscous flows, respectively. P, t], 0, are the pressure, 
which is calculated from the equation of state coiTesponding to 
a polytropic or to an ideal gas, the dynamical viscosity which is 
assumed to be identical to the shear viscosity, and (= V^u'^), 
which measures the divergence or convergence of the fluid world 
lines, respectively. W"' = u^'u" + g'^^ is the spatial projection ten- 
sor, whereas o" corresponds to the symmetric spatial shear ten- 
sor: 



For the X-ray burst calculations, the general relativistic Navier- 
Stokes equations are solved using the Boyer-Lindqueit coordi- 
nates in the background of a slowly rotating NS, with the fol- 
lowing metric elements: 
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V=g = p2 COS 6* = a Vt. 



In this formulation, the parameter "Qns" denotes the spin of the 
neutron star, which is taken to be much smaller than the break-up 
frequency. jS'^ is the frame-dragging frequency associated with 
the rotation of the NS: 



P^ = Qfd = [^]J*(^) 



2GM, 1^„2 ^ 
[^](^)Rns^^ns, 



(3) 



where J, = I»Qns, and I, - |MnsR^s is the moment of inertia 
of the NS. The parameters: c, Mns, G, rg(= ^^), a denote the 
speed of light, mass of the NS, the gravitational constant, the 
gravitational radius and the lapse function, respectively. In writ- 
ing these expressions, we made use of the coordinate transfor- 
mation = n/l - 0, where we use the latitude 9 instead of the 
polar distance angle 6; hence the appearance of "cos" instead of 
"sin" in the metric terms. 

The set of general relativistic Navier-Stokes equations in 3D axi- 
symmetry can be written as the residual vector equation: 



R = 0. 

The components of this vector read as follows: 

1 . The continuity equation 

R, = ^ + Ll.eD = 
ot 

2. The radial momentum equation 

R2 - ^ + LheM, - f,. - LllgM, = 

3. The vertical momentum equation 

R3 = ^ + LhgMg -fg- L2%Mg = 

4. The angular momentum equation 



(4) 



r/^" = V<,-u''h«"'' + V<;.u''h«"''. 



R4 = 



dt 



LlrflM^-f^-L2^„M. = 



(5) 
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5. The internal energy equation 



— + LheS." + Ll,flu'] = 0, 

ot at 



(9) 



decelerate, rather than accelerate, the convergence of the numer- 
ical solution procedure. In 3D axi-symmetry, few of these terms 
can be simplified or even neglected without violating the physi- 
cal consistency of the numerical scheme. 



where Ll,-g are first order advection operators that have the form: 

fr,e,ip are force terms that include pressure gradients, centrifu- 
gal and gravitational forces acting along the radial, horizon- 
tal and azimuthal directions, respectively. D(= pu') is the 
modified relativistic mass density. are the four-momenta: 

(Mt, Mr, Me, M^) = D(ut,Ur,Ufl,u^), where D = Dh, and u' is 
the time-like velocity, V - u^/u' is the transport velocity and 
"h" denotes the enthalpy. L2^g are the spatial projections of the 
viscous stress energy tensor T^^^, (see Eq. [1) in the respective 
direction. These are obtained from the projection of the viscous 
tensor along the vector normal to the hyperspace, i.e., constant 
in time: 



'^'^re - ^t^^ Vis ~ "i^^ Vis 



^ /ii-* Vis' 



Temperature 




Fig. 1. The starting distribution of the temperature in the domain of 
calculation. Low (high) values correspond to blue (yellow) color The 
red color corresponds to the hot bubble. 



^ corresponds to the spatial divergence of a 



where ^ = {r,Q,ip\. V 

tensor taken in the Boyer-Lindquist coordinates and F^^ are the 
Christoffel's symbols of the second kind 



For completeness, we re-write the forms of these second order 
viscous operators explicitly as follows: 
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To outline our simplification strategy, we first mention the fol- 
lowing relevant issues: 

1 . In most astrophysical fluid flows the molecular viscosity is 
too small to be relevant on observationally reasonable time 
scales. Thus, in the absence of other forms of viscosity, the 
above operators can be safely neglected. 

Moreover, these operators must vanish asymptotically when- 
ever the fluid velocity approaches the speed of light. 

2. Turbulent viscosity is more common in modeling astrophys- 
ical fluid dynamics. In rotating astrophysical flows, turbulent 
viscosity is a fundamentally important mechanism for angu- 
lar momentum transport. Therefore, the viscous operators of 
the angular momentum equation are important and should 
converge to the usual Newtonian form whenever the velocity 
becomes sub-relativistic. 

3. The viscous operators appearing in the radial and vertical 
momentum equations act, in general, to diffuse and smooth 
strong velocity-gradients. The mixed-derivative appearing 
in these two equations act to mainly enhance the viscous- 
coupling between the velocity components. 

Therefore, our simplification strategy relies on considering only 
second order, mixed-free and Laplace-like operators. In numeri- 
cal matrix algebra, such operators generally enhance the diago- 
nal dominance of the coefficient matrix and stabilize its inversion 
procedure. 

Specifically, using our time-implicit formulation, the following 



operators have been included into the numerical solver: 
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where is the contravariant form of the metric tensor 



2. 1 . Simplifying considerations 

Most of the above-described collection of viscous terms con- 
tains highly non-linear, first and second order operators, some 
of which are difficult to handle numerically or unnecessarily en- 
large the band-width of the coefficient matrix, while others may 



L2:, 



{ri [(f) (u'f 

■ W (w) - -3^0 ■ u') (ueu« + 1) + (§) (uV . 



1 



ae ■ 
, au„ 



.1 



= V, ■ 77 [(^ - i(g'^'%l + g'^^^f^)u^) (u'y + g-) 



flr 



-(i(g-^) u%u' - (i(g"%) u'^(u^u'^ + 1)] 



)]) 



(14) 



1 „90'>%' 



') 



(15) 



89 > ^ 



89 

5g''%^) u'^(u^u^ + 1)] 



4 



A. Hujeirat et al.: Low-Mach number flows 



T, t = 0.25 s 



1.4 - 
1.2 

o 

i 1.0 

T 
0) 

E 0.8 
S 0.6 

C 

o 

0.2 

0.0 



Grid distribution 




c 



0.0 0.2 0.4 0.6 0.8 i .u 
Radius 



Density, t = 0.5 s 



0.25 L 

0.95 1.00 1.05 1.10 1.15 1.2( 
Radius 

T, t = 0.5 s 




n 25 r , , , ,\ I , A , , i\, , , ,\ I , , A , I , \ , , 0.25 [ 

0.95 1.00 1.05 1.10 1.15 1.20 0.95 1.00 1.05 1.10 1.15 1 .2( 



Density, V, t = 0.75 s 



T, t = 0.75 s 





. J . .'. . I -I 0.25 L 

0.95 1.00 1.05 1.10 1.15 1.20 0.95 1.00 1.05 1.10 1.15 1.2( 



Density, t = 1 .0 s 




t = 1.0 s 



0.95 1.00 1.05 1.10 1.15 1.20 0.95 1.00 1.05 1.10 1.15 1.2( 



Fig. 2. a rising bubble in the white dwarf atmosphere. The domain of calculation is covered by 100 x 170 non-uniformly distributed finite volume 
cells (top left panel). The colored images are snapshots of the temperature of the rising bubble after 0.25 , 0.5. 0.75 and 1.0 second after the flash, 
low-to-high values of temperature correspond to blue-to-red. On the left panel, snapshot of the density distribution at 0.5, 0.75 and 1.0 seconds in 
black-white are shown. 



It can be easily verified that in the non-relativistic regime the where rj — pv and v denotes the kinematic viscosity, 
radial component of the diffusion operator L2)°^ reduces to the 
classical Newtonian form: 

= (16) 
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Fig. 3. The distribution of the angular velocity in the bursting region for four different viscosity parameters amr shortly after burst events (« 0.1ms). 
In the left panel 30 uniformly distributed isolines of the angular frequency fl are shown. On the right panel we display the radial profiles of £1 
across the bursting region. The soUd line in "(b2) corresponds to a relaxed fl-proflle after 10 ms. 



2.2. Viscosity prescription 



Therefore, we adopt the turbulent viscosity prescription: 

Vtur =< Vtur >< Aur >~ ^turVs X Ori/^NS, (17) 



Similar to classical accretion disks, we assume that molecular where < Vtur >, < ^tur > correspond to mean values of velocity 

viscosity is too small to have a significant effect on the angular and length scale of eddies in a turbulent medium, respectively, 
momentum distribution on a short time scale such as the ther- 
monuclear one. 
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These are set to be respectively smaller than the sound speed Vs 
and smaller than the radius of the NS. Thus, fftur, ai are con- 
stants that are set to be smaller than unity. In the present paper, 
all model calculations assume < ^tur 0.1/?ns, i-C-, ^2 = 0.1, 
whereas the parameter atur may differ from one model calcula- 
tion to another. 

3. Heat bubble calculations 

3.1. Numerical solution method 

The set of hydrodynamical equations are solved using a pre- 
conditioned defect-correction iteration procedure. The matrix 
equation to be solved in each iteration is: 

A/j = d, (18) 

where A is a preconditioner, i.e., a coefficient matrix that is sim- 
ilar to the Jacobian J, but which is much easier to invert. J is 
obtained by calculating the entries resulting from dR/dq, where 
R denotes the vector of equations (Eq. 4) and q the vector of 
variables. 

In this formulation jj = q'^^-q' corresponds to the correction be- 
tween two successive iterations and d is the defect (for a detailed 
description of the method see Hujeirat, 2005; Hujeirat et al., 
I2008h . 

We mention that the solver employed here relies on the con- 
servative formulation of the hydrodynamical equations, using 
the finite volume formulation. For strongly time-dependent sim- 
ulations, an advection scheme of thkd o rder spat i al and sec- 
ond order temporal accuracies is used (iHuieiratl l2005l) . As 
pre-conditioner we use the approximate factorization method 
(AFM), which is proven to be most appropriate for modeling 
weak ly and strongly low Mach number flows (IHuieirat et al.L 
I2007h . 

Rad. Velocity (6 = 0.33), a = 0.1, T = 3e 




1.05 1.10 1.15 1.30 1.25 1.30 
Radius 

Fig. 4. The profile of the radial velocity across the heat bubble shortly 
after the burst (0.3 ms), using o-tu, = 0.1. Obviously, the outward- 
oriented fluid motion is much stronger than the inwards one. 



3.2. Heat bubble propagation in the atmosphere of a 
non-rotating white dwarf 

Rising bubbles in stellar environments have been extensively 
studied by Almgren at al. (2006, see also the references therein). 
To test the capability of our solver at capturing the propagation 
of strongly time-dependent heat bubbles in strong gravitational 
fields, such as X-ray bursts on neutron stars, we adopt the same 



setup of the heat b ubble problem as described in Sec. 4.3 of 
l Almgren et all 12006). 

The domain of calculation is restricted to the first quadrant: 

D = [Ri^ < r < «ou,] X [0 < 6» < nil] 
= [1 <r< 1.35]x[0<6><7r/2], 

where length scales are measured in units of the radius of the 
star's core. 

The domain D is divided into non-uniformly distributed finite 
volume cells: 100 in the radial and 170 in the horizontal direc- 
tion, where the minimum grid spacings is set to coincide with 
the center of the initial heat bubble. 

The initial conditions are constructed by carrying out 3D axi- 
symmetric time-dependent calculation to obtain the hydrostatic 
equilibrium. A small region is then chosen (Fig. [1]), where the 
matter is replaced by a thinner but much hotter plasma, while 
keeping it in pressure equilibrium with the surrounding media. 
In Fig. |2] we show several snapshots of the rising bubble. 
Obviously, these results agree quite well with the unsplit and the 
low Mach number schemes presented bv Almgren et al.. (.2006) . 

3.3. Rotating heat bubbles in deep gravitational fields of 
neutron stars 

Similar to the Sec. (3.2), we apply our general relativistic solver 
to model the rise of a rotating bubble starting from below the 
atmosphere of a rotating neutron star. The radius of the neutron 
star's core /?ns is set to be equal to the following radius, 

7?NS = 3rg(l-f ^l-02^), (19) 

which is smaller than the classical last stable orbit of a 
Schwarzschild black hole. The inner boundary of the domain 
of calculation is taken to be the radius of the core, whereas the 
outer boundary is located at 7?out - 1.35/?ns- The core is set to 
rotate at 0.4 the break-up velocity, whereas Q at the outer radius 
is 10 times smaller. The matter in the domain is set to adjust its 
rotation to the boundaries through viscous interaction. 
The initial distribution of the variables are obtained by solving 
the general relativistic Navier-Stokes equations to obtain station- 
ary and differentially rotating flow configurations. As in the pre- 
vious section, the heat bubble is injected into the domain through 
replacing the medium at a certain location by a hot and tenuous 
plasma, while keeping it in pressure and rotational equilibrium 
with the surrounding media. 

Unlike the calculations in the pervious section, the purpose of 
the present calculations is to unveil the response of the surround- 
ing region to violent events associated with a dramatic change in 
the distribution of the angular momentum. 
Therefore we run several numerical calculations with 
fftur - 1.0, 0.1, 0.01, 0.001. The initial stationary con- 
figurations have been obtained using the corresponding value of 
atur- The results are shown in Fig. [3] and can be summarized as 
follows: 

1. All model calculations show a pronounced deficiency of an- 
gular momentum at the central part of the burst, accompa- 
nied with a significant increase at the boundary of the bub- 
ble. Thus, the burst leads to the formation of a dynamically 
unstable flow-configuration: a shell of slow rotating matter is 
bounded both from below and from above by relatively fast 
rotating matter. 
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Fig. 5. 25 isolines of the angular velocity (left) and of the Lorentz factor (right panel) overlayed on their corresponding colored images after 0.5 
ms are shown, using o-tu, = 0.01. Blue-to-red color correspond to low-to-high values of CI an T. 



2. The viscous-induced fronts of angular momentum are found 
to propagate outwards at much faster speed than inwards. 
This implies that the matter in the deeper layers adapts its 
conditions to the inner boundary much faster than at the outer 
boundary. This is a consequence of the adopted prescrip- 
tion of the viscosity r], which is more effective in hotter and 
denser regions of the plasma. Also, the outwards-oriented 
velocity is obviously larger than the inward one (Fig.|4]i. 

3. The difference between the rotational velocity at the center 
of the bubble and that at its boundaries becomes more sig- 
nificant, the smaller the atur-parameter is chosen. 

In the fftui- = 1 .0 case (see Fig.[3](al)), the rise time of the bubble 
is found to be extremely long compared to those obtained with 
smaller atur- 

Here the viscous pressure in the radial direction P^^,^ ~ '^tui-|7 
has an opposite sign and acts to reduce the effects of the thermal 
pressure Pth- In the extreme case, when f th + ^"¥180 — > 0, the 
effective sound speed, V^-^^ ~ 5(f the + Pvisc)/Sp — > 0, hence the 
propagation time rpro ~ R/V^-^^ — > oo. 

On the other hand, Tpro becomes of the order of one second, when 
using atur ~ <5(10 '), which fits well into the observed duration 
of burst events on NSs. 

In the low atur cases, (see Fig. [3]bl,cl,dl and Fig.|5]l, the sig- 
nificant increase of the rotational velocity at the outer O-fronts 
is obvious, but unreasonably large. In the Fig.|2](dl), the matter 
at the outer front is found to be gravitationally unbound to the 
central NS, giving rise to strong outflows. 
However, as outflows during X-ray bursts can be excluded on 
observational grounds, we conclude that atur must acquire much 
larger values, and that specifically atur ~ (9(10"'). 
To study the viscous-reaction of the matter in the adjusting lay- 
ers to the sudden increase of rotation induced by a burst, we have 
run separate calculations in the following manner A solution for 
the hydrodynamical equations including rotation has been hy- 
drodynamically calculated. As a second step, we modified the 
Q-profile by including a Gaussian perturbation of the form de- 
picted in Fig. ^tq. We then followed the time evolution of this 



profile on the viscous time scale. The profiles ri,T2,Tio corre- 
spond to Tvisc/lO, -Tvisc/S, ^visc ■ 

These calculations show that the effect of turbulent viscosity is 
to mainly transport angular momentum outwards. As a conse- 
quence, the deficiency in the rotational support in the deep lay- 
ers enhances the compression of matter and give rise to an addi- 
tional burst in the neighboring shell. This chain of reactions may 
run away to spread over the whole surface of the NS on the vis- 
cous time scale, which is of the order of one second, assuming 
atur -0.1. 

When the outer layers cool, the turbulent viscosity decreases and 
the corresponding viscous time scale increases as Tvisc ~ ^HT . 
This implies that after the burst, the time scale needed for the 
outer layers to adjust their rotation to the bulk of the star might 
lengthen by an additional order of magnitude. As a consequence, 
the observed spin up of NSs in their cooling tail is a manifesta- 
tion of the increased rotational velocity of the outer layers caused 
by the burst events, but which, eventually, should decrease at 
later times. 

4. Summary & Conclusions 

In this paper we have presented the set of general relativistic 
Navier-Stokes equations in 3D axi-symmetry, using the Boyer- 
Lindquist coordinates in the background of a slowly rotating 
neutron star. We have shown that, for modeling relativistic vis- 
cous flows, the collection of viscous operators can be reduced 
considerably into a smaller class, which consists of the domi- 
nant second order operators, and which subsequently enhances 
convergence of the implicit solution procedure. 
The set of equations are solved numerically using an implicit 
solution procedure, which is based on a pre-conditioned defect- 
correction iterative method. Similar to Taylor-flows between 
concentric spheres, we use the "Approximate Factorization 
Method" as a pre-conditioner, which has superior converging 
properties over other non-symmetric methods, such as the black- 
white line Gauss-Seidel method. In the non-rotating case, we 
have shown that the solver is capable of reproducing the time- 
evolution of heat bubbles during their rise to the surface of 
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Fig. 6. Successive snapshots of the omega-profiles in the outer layers of 
a neutron star. Shell I shows the burst-induced propagation of Omega- 
fronts into the neighboring shells. In shell II, the viscous-induced reac- 
tion to the Omega-fronts is shown. This schematic picture shows that 
viscosity transport the excess of angular momentum preferably in the 
outward direction. The lower figure, which is obtained using hydrody- 
namical calculations, clearly confirms this behavior. 

the wh ite dwarfs, as reported by I Almgren et al.l (l2006l) : lNonakal 
(I2OO8I) . accurately. 

In the rotating case, it has been shown that viscous-generated 
fronts inside heat bubbles propagate into the surrounding quite 
rapidly. The effect of these front is mainly to transport angular 
momentum to the outer layers, leaving the matter in the deeper 
layer with less rotational support, hence more compressed. 
The increased rotational velocity of the outer layers may be con- 
nected to the observed spin up of neutron stars during X-ray 
bursts. However, this increase will eventually be followed by a 
spin down on later times when the outer layers cooled down to 
their pre-burst thermal state. 

Our numerical results show that a viscosity parameter of the or- 
der of atm- ~ 0. 1 is the most suitable value for fitting observa- 
tions of NSs during X-ray bursts. A much larger value yields 
propagational time that is much larger than one second, whereas 
smaller values yield unstable shell configurations and gives rise 
to gravitationally unbound outflows. 

In addition, a possible mechanism that may underly the 
rapid spreading of burning ignition fronts has been presented. 
Accordingly, the viscous-generated fronts inside a heat bubble 
may transport angular momentum into the horizontally adjusting 
layers. The viscosity then acts to decouple matter from angular 
momentum, subsequently enhancing compression in the deeper 
layers and giving rise to thermonuclear runaway. This chain of 
reactions may run away to spread over the whole surface of the 
NSs on the viscous time scale Tyis ~ 1 s . 

Concerning frame dragging, while included in the present sim- 
ulations, its effects are still to be quantified. Also, the effects of 



magnetic fields and thermal diffusion is the subject of an ongo- 
ing work. 
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